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Abstract. We present a calculation of the next-to-leading order cross section for the inclusive 
hadroproduction of D and B mesons as a function of the transverse momentum and the rapidity 
in a massive variable flavor number scheme. We compare our numerical results with recent data 
from the CDF Collaboration at the Fermilab Tevatron for the production of Z)°, D* + , D + , and D+ 
mesons at center-of-mass energy \/S = 1 .96 TeV and find reasonably good agreement with the 
measured cross sections. 
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Various approaches for next-to-leading order (NLO) calculations in perturbative QCD 
have been applied to one-particle inclusive hadroproduction of D or B mesons. For 
definiteness, we shall consider here D mesons. However, all results can easily be carried 
over to any other heavy-flavored hadron. 

A basic approach is the fixed flavor number scheme (FFNS) (H), in which the number 
of active flavors in the initial state is fixed to rif = 3 and the charm quark appears only 
in the final state. The charm mass m is explicitly taken into account together with the 
transverse momentum pj of the observed meson. In this scheme the charm mass acts 
as a cutoff for the initial- and final-state collinear singularities and collinear logarithms 
ln(/?y/m 2 ) are kept in the hard scattering cross sections. However, for pj 3> m, these 
logarithms become large and spoil the convergence of the perturbation series. 

Therefore, in the regime pj ^> m, it is more appropriate to treat charm quarks like 
massless partons and to absorb the collinear logarithms into scale dependent parton dis- 
tribution functions (PDFs) and fragmentation functions (FFs). As is well-known, by this 
procedure the large logarithms ln(p^/m 2 ) are summed via the DGLAP evolution equa- 
tions and the hard scattering cross sections are finite (infrared safe) in the limit m — > 0. If 
the power-like charm mass terms 0(m 2 / pj) are neglected this is just the conventional 
parton model or zero-mass variable flavor number scheme (ZM-VFNS). Usually, in the 
ZM-VFNS the charm mass is neglected from the beginning and the collinear singulari- 
ties appear in dimensional regularization as poles in e where d = 4 — 2e is the number 
of space-time dimensions. Conventionally, these poles are removed in the modified- 
minimal- subtraction (MS) scheme. If, on the other hand, the collinear singularities have 
been regularized with help of a mass m it is necessary also to subtract finite terms along 
with the collinear logarithms lnm 2 in order to recover the hard scattering cross sections 
in the MS scheme. 
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FIGURE 1. QCD predictions for one-particle inclusive production of charmed mesons, X c = 
D° , D* + , D + , D+ , at the Tevatron Run II. In each case, the results are shown for the average of the observed 
meson with its antiparticle (X c +X c )/2. The solid lines have been obtained with /Zr = jj,p = jj, F = mj. The 
upper and lower dashed curves represent the maximum and minimum cross sections found by varying /Xr, 
/If, and Hp independently within a factor of 2 up and down relative to the central values while keeping 
their ratios 0.5 < /If / Hr, / Mr, / Mf — 2- CDF data |9] are shown for comparison. 



On top of these two basic approaches, schemes have been devised which combine 
the two features, non-zero charm mass and resummation of ln(p^/m 2 )-terms. One such 
scheme, which has been applied already to inclusive charmed meson production for 
the Tevatron experiment is the so-called fixed-order next-to-leading-logarithm (FONLL) 
scheme. This scheme smoothly interpolates between the traditional cross section in the 
FFNS and a suitably modified cross section in the ZM-VFNS approach with perturbative 
FFs with the help of a pj dependent weight function ylyD. In both non-zero-charm- 
mass approaches, FFNS and FONLL, the theoretically calculated cross sections are 
convoluted with a scale-independent non-perturbative FF extracted from e + e~ data 
describing the transition from the produced charm quark to the observed D meson. 

Recently, a general mass variable flavor number scheme (GM-VFNS) has been 
worked out by us 0, El 0] which is closely related to the ZM-VFNS, but keeps all 
m 2 / p\ terms in the hard- scattering cross sections in order to achieve better accuracy 
in the intermediate region pj > m. The massive hard scattering cross sections have 
been constructed in a way that the conventional hard scattering cross sections in the 
MS scheme are recovered in the limit pj — > °° (or m — > 0). The requirement to adjust 
the massive theory to the ZM-VFNS with MS subtraction is necessary since all com- 
monly used PDFs and FFs for heavy flavors are defined in this particular scheme. In this 
sense this subtraction scheme is a consistent extension of the conventional ZM-VFNS 
for including charm-quark mass effects. It should be noted that our implementation of a 
GM-VFNS is similar to the ACOT scheme which has been extended to 1 -particle inclu- 
sive production of B mesons a few years ago M\- There are small differences concerning 
the collinear subtraction terms o]. Further, in L8], the resummation of the final state 
collinear logarithms has been performed only to leading logarithmic accuracy. 

To calculate the cross section d 2 o /dprdy for the reactions p + p — ► D + X, FFs 
are needed which describe the fragmentation of the charm quarks, the light quarks, 
and the gluon into the observed D mesons. Fragmentation functions for the D* meson 
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FIGURE 2. Ratios of experimental results and our central theoretical predictions (solid lines in Fig. [Q. 
In addition, the theoretical uncertainty bands are shown, obtained as the ratio of the upper (lower) QCD 
prediction and the central curve. 



have been extracted at leading and next-to-leading order already some time ago iflpll . 
using experimental data from the OPAL [[Till and ALEPH lfl2ll collaborations at LEP1. 
Recently, using the same procedure as in fllOQ , also the FFs for Xr — D° ,D + ,Df : A r have 
been determined 0] using OPAL data for <?+e X c + X Q. In Refs. J3 El the 
fits for the X c FFs have been performed using a starting scale /io = 2m for the gluon and 
the w, d, s and c quarks and their antiquarks, while /io = 2m/, (m^ = 5 GeV) was chosen 
for the FFs of the bottom quark and antiquark. The FFs of the gluon and the first three 
flavors were assumed to be zero at this /io- At larger scale /i, these FFs are generated 
through the usual DGLAP evolution. Since the effect of the gluon FF is important at 
Tevatron energies as was found for D* production in [0] we decided to repeat the fits 
for the X c FFs with the lower starting scales /io = m and /io = m/„ respectively. This 
changes the FFs of the c quark only marginally but has a sizable effect on the gluon FF. 
The details of these new FFs will be presented elsewhere Q. 

Next we show our predictions for the cross sections do/dpr for D°,D* + 1 D + and 
Df production obtained in the GM-VFNS. For a comparison with the ZM-VFNS we 
refer to Ref. lflr3l . The partonic cross sections are convoluted with the (anti-)proton 
PDFs and the FFs for c -> X c , u, d : s -> X c and g ->X> We use CTEQ6M PDFs JL1 
and the FF sets for D°, D* + , D + and Df from H 1 311 . 1 Results are shown for the 
average of the observed X c mesons with their antiparticles. We consider do /dpr at 
v 7 ^ = 196 TeV as a function of pj with y integrated over the range — 1.0 < y < 1.0. 
For the charm mass we take m = 1.5 GeV and evaluate aj"^(/i#) with nf = 4 and 

scale parameter A^ = 328 MeV, corresponding to a^imz) = 0.1181. The results 
are presented in Figs. [T] and El The solid lines correspond to the central scale choice 
jlR = l±F = l±' F = mj = (pj -r-m 2 ) 1 / 2 , where /i R is the renormalization, /if the initial- 
state and n' F the final-state factorization scale, respectively. To investigate the scale 



It should be noted that the results presented at the DIS05 have been obtained with the FFs from 1 10 Hl3ll . 



variation of our predictions, we independently vary the renormalization and factorization 
scales by a factor of two: 0.5 < Uu/mr, Uv/mr, U^/mr < 2 while keeping their ratios 
0.5 < Hf/iir,ii' f /iir,ixf/ii' f < 2 HI 611 . Our theoretical results are compared with the 
experimental data from CDF [9]. As can be seen, the data are in good agreement with 
the upper curve of the uncertainty band whereas they are a factor of about 1 .5(1 .2) above 
our central prediction at low(high) pj. 

Residual sources of theoretical uncertainty include the variations of the charm mass 
and the assumed PDF and FF sets. A variation of the value of the charm-quark mass 
does not contribute much to the theoretical uncertainty. Also the use of other up-to-date 
NLO proton PDF sets produces only minor differences. Concerning the choice of the 
NLO FF sets we obtain results reduced by a factor of 1.2-1 .3 when we use the NLO sets 
obtained by fitting with the initial scale choice /lo = 2m, 2m^. 

In conclusion, we have presented a NLOperturbative QCD calculation of D meson 
production at the Tevatron in a GM-VFNS |4, 5] which provides the best description of 
these experimental results obtained so far. It completes earlier work in this scheme on 
D meson production in 77 and yp collisions lfl8ll . This approach will be applied next to 
B meson production at the Tevatron. Furthermore, it is planned to extend this scheme to 
heavy meson production in deep inelastic scattering. 
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